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We present a determination of a new class of three-loop Feynman diagrams describing heavy-to- 
", light transitions. We apply it to find the O^a^j corrections to the top quark decay t — > bW and 

f~>) ■ to the distribution of lepton invariant mass in the semileptonic b quark decay b — > ulv. We also 

' confirm the previously determined total rate of that process as well as the Oyot 2 ) corrections to the 

. muon lifetime. 
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The determination of higher order corrections in perturbative quantum field theory is notoriously difficult, and 
with the general tendency towards precision measurements in particle physics, each newly- won class of perturbative 
integrals expands the possibilities for phenomcnological analyses. For instance, quantum corrections to decays of 
I i neutral particles, such as a virtual photon or a Z boson into hadrons, are known to sixth order in perturbation theory, 
> ; <D(a 3 s ), and even some O(cvf) effects have been studied. Those results have been very useful in determining a variety 
t-H , of Standard Model parameters such as the Z boson properties, the strong coupling constant, and the running of the 
CN ■ electromagnetic coupling constant 0. 

Much less is known about radiative corrections to processes with a charged particle in the initial state. Only 
relatively recently have first results been obtained in fourth order perturbation theory, O(o^) and C(a 2 ), primarily 
for total decay rates. The technical challenge in such calculations is the presence of massive propagators. For example, 
consider the muon decay. Since the muon is charged, it can emit photons, and the resulting amplitudes will involve 
propagators of a virtual muon. Its mass sets the energy scale of the process and cannot be treated as a small parameter. 
The presence of massive propagators is an obstacle in evaluating the multi-loop diagrams required by precise 
II ' measurements of heavy quark and lepton decays. So far, genuine ©(o^) corrections to heavy quark decays are known 
Oh, only for semileptonic processes, Q — > qlv(gg), and only for some kinematic cases. One approach that has been 
successful consists in expanding Feynman diagrams around the zero recoil limit: when the quark q remains at rest 
with respect to Q. The kinematics of semileptonic decays can be represented by a triangle, since the invariant mass 
■ of the leptons together with the mass of the final state quark q cannot exceed the mass of the decaying quark. This is 
k^j ' depicted in Fig.^ The diagonal boundary corresponds to the zero recoil limit, in which the O(cVg) effects are known 
13- Also shown are the starting points of previously studied expansions. Those results have helped improve the 
d ' knowledge of the b quark lifetime and the determination of the CKM matrix element V c b- 

However, the expansion around zero recoil converges slowly near the origin in Fig. ^ that is, when both the quark 
q and the lepton pair are light; in this case computations become prohibitively expensive. Two other approaches have 
been used in such cases. First, for the phenomenologically important decays fi — > evv and b — > ulv, the total lifetimes 
have been determined in 0, |g by analytically calculating imaginary parts of four-loop diagrams. That heroic effort 
is difficult to extend, for example, to differential distributions. The other approach consists in expanding diagrams in 
an artificial parameter, for example the ratio of muon masses outside and inside loops, and using Pade approximants 
to sum the expansion for the physical value of this parameter. This approach was used to check the muon and b — ► u 
results 0, 0- It yields results scattered around the exact values which are often sufficient for applications. The 
drawback of this method is that it is very difficult to estimate the errors reliably. 

The purpose of this study is to extend the method of expansions beyond the zero recoil limit. We start directly at 
the origin of Fig. ^ which corresponds to the kinematics of a top quark decay into a massless b quark and a massless 
W boson. We then treat the W mass as a small perturbation and compute several terms of the resulting expansion. 
We stop when we can smoothly match to the previously obtained expansion around the case of the W boson equally 
heavy as the decaying quark |!J]. The physical top quark decay corresponds to a specific value of the W mass, but in 
the more general decay Q — > q + leptons, the term "W mass" refers to the invariant mass of the lepton pair, and it is 
in this context that we can relate t — > bW to b — > ulv and u — > evV. 

This is the first time that exact results are available in this limit and we can now address a number of interesting 
problems. We obtain an accurate value of the O{oq) correction to the top quark lifetime. The combination of 
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FIG. 1: Kinematical boundaries of the semileptonic decays Q — > q + leptons. The solid arrow shows the expansion presented 
in this paper. Previously known expansions are indicated with dotted arrows. 

our results with the expansion around the heavy W case allows us to give a complete description of the differential 
distribution of b — > ulv decay in the invariant mass of leptons, and thus improves the theoretical description of this 
decay, important for the determination of V u b- We also check the muon and b — > u lifetime corrections with a relative 
error of about 2 x 10~ 4 . In the future, the same method can be employed to improve perturbative corrections to 
mixing processes such us Bd <-* Bd and B s <-» B s . 

In Fig.[21we show three examples of the diagrams that we have to consider in order to calculate t — > bW at 0(a 2 ^). 
We use the optical theorem to connect the imaginary parts of such diagrams with contributions to the decay. Note 
that we customarily speak about two-loop corrections when what we actually need to compute are the imaginary 
parts of three-loop diagrams. The various cuts correspond to two-loop virtual corrections or emissions of one or two 
real quanta. 



FIG. 2: Examples of diagrams whose cuts contribute to the O(af) decay rate for t — > bW: (a) light or heavy quarks; (b) 
abelian; (c) non-abelian. 

With mj, = 0, there are two scales in the problem: m t and raw- We define an expansion parameter to — jm\ so 
that the two scales can be expressed as hard and soft (0(1) and O(yfuj), using mt as the unit of energy). Contributions 
arising from these two scales are identified using asymptotic expansions so that we must consider two regions. In the 
first region, all the loop momenta are hard and the W propagator can be expanded as a series, in powers of lo, of 
massless propagators. In the second region, the gluon momenta are hard but the loop momentum flowing through 
the W is soft. In this region, the diagrams factor into a product of a two-loop self-energy type integral and a one-loop 
vacuum bubble integral with a scale of m\y. The leading contribution from this second region is C(w 2 ), and the 
interplay between the two regions gives rise to terms with a large logarithm lnw. 

All scalar integrals arising in the problem can be expressed in terms of 9 basic topologies. We use differential- 
algebraic identities to reduce all loop integrals in both regions to a combination of 24 master integrals. The resulting 
large linear systems can be solved in a few ways. In the traditional method [Tcj . one inspects the structure of the 
identities and rearranges them manually into the form of recurrence relations for an efficient iterative solution of the 
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system. This "by inspection" method has proven to be very successful in numerous applications (e.g., 0,0, 0, |n|) but 
it requires much human work to implement. Conversely, a straightforward solution of the linear system is much more 
expensive computationally and was first achieved only recently ^ n our calculation 0] we used the traditional 
approach (programmed in FORM [T^jl as well as a modified version of the new algorithm for which we implemented a 
dedicated computer algebra system. In both cases we independently obtained identical results which serve as a check 
of correctness but also enable us to compare these two methods. Details of the implementation of both methods and 
the evaluation of master integrals will be presented in a forthcoming technical paper. 
The final result for the top quark decay width can be written as 



T(t -> bW) = r 
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Throughout this paper, we use a s = a 
0(a s ) coefficients are already known analytically |15| 
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' IS ([i), where /j, is the pole mass of the decaying quark. The tree- level and 
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The C>(a 2 ) result can be subdivided into four gauge-invariant pieces, 

X 2 = Cf (T r N l X l + T R N H X H + CfXa + CaXna) , 



(2) 
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where Cf ■ 
light (m q = 
series to at 



= 4/3, Ca = 3, and Tr = 1/2 are the usual SU(3) color factors and Nl and Nh denote the number of 
0) and heavy (m q = m t ) quark species. For the coefficients Xl, Xh, Xa, and Xna, we have obtained a 
least u! 5 , of which the leading terms are 
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The leading term, O(ui ), of these results can be compared with the numerical estimates obtained with the zero 
recoil expansions in Eq. (14) of [3; all of our results agree within their error estimations. Our result can also be 
compared with a numerical study of the top decay rate obtained by means of Pade approximations up to O(to 2 ) 0. 
In many cases we find agreement. However, there are also instances where the numerical estimates in differ from 
our analytic expressions JSJl by a few error bar lengths, illustrating limitations of the Pade approximation in this 
problem. For example, the coefficient of the 0(ui) term of the nonabelian part Xna of Eq. 101 is 3.3398 whereas the 
value cited in [j| reads 3.356(3), corresponding to a 5<r discrepancy. Similarly in Xa, the O{uo) term is off by 3<r. 

With a sufficient number of terms, the present expansion can be smoothly matched with the one around the uj — 1 
limit studied previously [9j in the context of semileptonic b quark decays. The result of such a matching procedure 
is depicted in the graphs in Fig. |21 Although strict matching of the two expansions in the entire interval < u> < 1 
would require a very large number of terms from each side, a wide overlap region arises even when only a few terms 
are taken into account. 

The most obvious application of the above result is the precise determination of second order QCD corrections to 
the top quark decay rate. An estimation of this effect is already known, both from numerical studies and from an 
extrapolation of the zero recoil limit. However, for the measured ratio of W and top masses, to ~ 0.213 [l7], the 
present expansion is the best way to calculate an accurate value of this contribution with a reliable error estimate. 
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Our expansion gives X% = —15.5(1) where the uncertainty is almost entirely due to the experimental uncertainty of 
m t . The theoretical error, which originates from taking a finite number of terms in our expansion, is 20 times smaller 
and can be still easily reduced if needed. Using a s (mt) — 0.11, we find that the two-loop correction decreases the tree 
level decay rate by about 2%, in agreement with earlier expectations. 

Our result also provides a check of the total lifetime calculations carried out for /i — > evu and b — > ulv decays. 
In these processes the expansion parameter oj corresponds to the invariant mass of leptons produced in the decay 
and our matching procedure allows us to obtain a differential width dT/dio valid in the full range of u with desired 
accuracy. The inclusive semileptonic decay rate b — ► ulv can be calculated by integrating over u> within the kinematical 
boundaries. Taking Ni — 4 and Nh = 1, we end up with du;X2{uj) — —10.644, which almost perfectly reproduces 
the —10.648 given in |6J. Analogously, the two-photon correction to the muon lifetime emerges from an integration 
of the abelian contribution Xa- We find J Q <1ujXa(u) = 1.7797, which is in excellent agreement with the exact result 
1.7794 @. 

To summarize, we have presented a new analytic 0{a 2 ^j result for the decay t — > bW in terms of a parameter 
u> = m^/m 2 and in the limit of raj, = 0, corresponding to the last remaining kinematic region in which the 0(a%\ 
heavy quark decay rates were not analytically known. This result has enabled us to confirm or modify slightly 
the corresponding results of previous numerical calculations. Our formulas are readily applicable to other physical 
processes such as muon decay and the semileptonic b quark decay b — > ulv. 

Our results depend on the imaginary parts of a novel class of three- loop integrals, which we have obtained using 
two independent paradigms for the solution of large systems of recurrence relations. To the best of our knowledge, 
this is the first time that both approaches have been used simultaneously to obtain a new result, and an objective 
analysis of the strengths and weaknesses of each approach will increase the efficiency of other large calculations in the 
future. This augurs well for the increasingly difficult physical problems that lie ahead. In particular, the top quark 
decay problem considered here has laid the foundation for 0(a 2 ^) perturbative calculations of mixing processes such 
us Bd <-* Bd and B s <-> B s . Since the recently found 0(a s ) effects are large and suffer from strong scale dependence, 
such improvement will help use those processes as a probe for new physics. 
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